Fermi surface dichotomy in systems with fluctuating order 
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We investigate the effect of a dynamical collective mode coupled with quasiparticles at specific 
wavevectors only. This coupling describes the incipient tendency to order and produces shadow 
spectral features at high energies, while leaving essentially untouched the low energy quasiparticles. 
This allows to interpret seemingly contradictory experiments on underdoped cuprates, where many 
converging evidences indicate the presence of charge (stripe or checkerboard) order, which remains 
instead elusive in the Fermi surface obtained from angle-resolved photoemission experiments. 
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I. INTRODUCTION 

Strongly correlated systems like the heavy fermions 
and the superconducting cuprates display excitations 
over a variety of energy scales. On short timescales (high 
energies) electrons are excited incoherently over energies 
ranging from the highest local (Hubbard U) repulsion to 
the magnetic superexchange interaction. On the other 
hand, over long timescales, i.e. at low energies around the 
Fermi level, the excitations can acquire a coherent char- 
acter typical of the long-lived Fermi-liquid quasiparticles 
(QPs). These energy scales appear very clearly both in 
theoretical^ and experimental 2 - studies of the one particle 
spectral function of strongly correlated systems. It is usu- 
ally assumed that the incoherent part has no momentum 
structure, an assumption which is emphasized by infinite 
dimensional studies where the self-energy is momentum 
independent and spatial informations on quasi long range 
correlations close to a phase transition are lost pi 

In this paper we discuss how this picture is modified in 
physical spatial dimensions. We want to address how the 
spectral function looks like when the system is close to 
an ordered phase. This issue is particularly important in 
the context of cuprates where it has been proposed that 
some kind of stripe-like order fluctuates in the metallic 
phase^i 4 ^ 

The scenario that we propose is based on the following 
qualitative argument: In physical dimensions a system 
may have long (but finite) ranged order parameter spa- 
tial correlations which are also long lived close to a quan- 
tum critical point. This defines a fluctuating frequency 
u>o above which the systems appears to be ordered. We 
argue that for energies larger than ujq with respect to the 
Fermi level the spectra should resemble the spectral func- 
tion of an ordered system. This spectral weight resides 
in what is usually called the incoherent part, which we 
argue, can have some important momentum structure. 
On the other hand at lower energies electrons average 
over the order parameter fluctuations and "sense" a dis- 
ordered system. In this limit we expect Fermi-liquid QPs 
with all their well known characteristics like a Luttinger 



Fermi surface (FS). 

To understand the momentum structure of the spectral 
function at energies higher than loq is important because 
if the incoherent part carries a memory of the close-by 
ordered phase it should be possible to analyze it to ob- 
tain informations on what is the underlying fluctuating 
order. Usually ordered systems are well described by 
mean-field thus one can obtain a first guess of how the 
incoherent part of the spectral function in the disordered 
phase should look like by performing a mean-field com- 
putation assuming long-range order. Comparison with 
experimental data in the absence of long range order can 
be useful to identify the fluctuating order parameter. 

To fix ideas consider as an example a moderately large 
U Hubbard system in a half-filled bipartite lattice in two 
dimensions at T = 0. In this case an antifcrromagnetic 
state is expected to be a competitive low-energy state. 
When the system is in the ordered phase the spectral 
function will be reasonably well described by a mean-field 
computation and will show two Hubbard bands separated 
in energy by mU with m the staggered magnetization. 
The bands will show some dispersion governed by the 
scattering of the electrons with the mean-field potential. 
Suppose that due to some frustrating effect long range 
magnetic order is lost while keeping well formed magnetic 
moments. We expect that beyond mean field if U is not 
too large (so that the disordered phase is metallic) a QP 
will appear at the center of the Hubbard bands with small 
spectral weight resembling the dynamical mean field the- 
ory (DMFT) pictured At high energies, however, elec- 
trons will sense a mean-field-like staggered potential for 
distances of the order of the correlation length £, which 
can be quite long, and therefore the system will keep sub- 
stantial memory of the mean-field like bands with their 
dispersion. Roughly we expect that the spectral func- 
tion will look like the superposition of a Fermi-liquid-likc 
spectral function, with a small weight z close to the Fermi 
level, plus a blurred mean-field-like spectral function in 
the presence of long-range order with a large weight 1 — z . 
This is at first sight similar to the DMFT picture but it 
differs in that in DMFT there are no magnetic correla- 
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tions surviving in the disordered phase and the incoher- 
ent part becomes momentum independent. We will show 
that in finite dimensions the incoherent part carries im- 
portant informations encoded in the momentum depen- 
dence which within a DMFT approach would require the 
cluster extensions developed more recently.— 
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FIG. 1: Schematic view of heavy-fermion system with a 
Kondo-like resonance arising at the Fermi energy E f from the 
mixing of a deep narrow / level (not shown) and the conduc- 
tion band (dashed line). Two QP bands arise E 1 ' 2 {k). The 
corresponding momentum distribution function rik is shown 
below with a true Fermi momentum kp P and a "fictitious" 
Fermi surface at k F . 



Another example can clarify the concept of an inco- 
herent part with a strong momentum dependence which 
carries physical information on the short range physics. 
Lets consider the more standard issue of large and small 
FSs in heavy fermions represented in Fig. [TJ In heavy 
fermions strongly correlated electrons in a narrow half- 
filled / level hybridize with electrons in a conduction 
band and give rise to a Kondo resonance at the Fermi 
level formed by coherent QP states. The width (and 
weight) of this QP band is usually quite small and sets 
the scale of the coherence energy in these systems. Now 
consider the momentum distribution function defined by: 



duM(k,u;)/M, 



(1) 



where f(u>) is the Fermi function and the spectral density 



S"(k,w) 



A(k, w) = -ImG(k, w) = - 

7T 7r(cj-S'(k,w)) 2 +E" 2 (k,w) 

is proportional to the imaginary part of the electron 
Green function with real (imaginary) part of the self- 
energy £' (S"). It is crucial to recognize that rtk involves 



all the excitation energies and its features might be dom- 
inated by the incoherent part of the spectrum if the QPs 
have a minor weight. Indeed strictly speaking the true 
FS at zero temperature is given by the small jump in the 
Fermi distribution function determining the Fermi mo- 
menta of the QPs at k^ P . This FS is large and satisfies 
the Luttinger theorem with a number of carriers includ- 
ing the electrons in the / level. This FS would naturally 
be determined by following the QP dispersion. On the 
other hand the shape of nk is substantially determined 
by the (incoherent) part of the spectral function, which 
has strong weight at energies corresponding to both the 
/ level and the conduction band. This latter gives rise to 
a rather sharp decrease of 7ik at a "fictitious" Fermi mo- 
mentum k c , corresponding to the FS that the electrons 
in the conduction band would have in the absence of mix- 
ing with the f-level. If the hybridization between the / 
level and the conduction band is turned off so that the 
QP weight z is driven to zero one reaches a situation in 
which the decrease at k c becomes a discontinuity and the 
small jump at k^ P disappears. It is clear that the sharp 
decrease of spectral weight at k c for finite hybridization 
has strong physical content for an observer who ignores 
the underlying model. Mutatis mutandis this shows that 
a computation in which the fluctuating order is artifi- 
cially frozen can give some hints on the distribution of 
the incoherent spectral weight in the less trivial case with 
fluctuations. 

The above ideas but with a more complicated order pa- 
rameter may explain perplexing data on cuprates. Sev- 
eral years ago underdoped La2- a; Sr 2 ;Cu04 (LSCO) com- 
pounds were examined and a dichotomy was found in 
the Fermi surface (FS) determined by two different treat- 
ments of the data. On the one hand the momentum de- 
pendence of the low-energy part of the energy distribu- 
tion curves was followed, thereby reconstructing the low- 
energy QP dispersion. In this way a large FS was found 
corresponding to the Fermi-liquid LDA band-structure 
and fulfilling the Luttinger requirement that the volume 
of the FS encircled the whole number of fermionic carri- 
ers n = 1 — x. On the other hand the FS was determined 
from the momentum distribution n^, obtained by inte- 
grating the spectral function over a broad energy window 
(~ 300 meV). Then the locus of momentum-space points 
where displays a sharp decrease, marked a FS formed 
by two nearly parallel (weakly modulated) lines along 
the k x direction and crossing two similar lines along the 
k y direction. This crossed FS would naturally arise in a 
system with one-dimensional stripes along the x and y 
directions. 

As a matter of fact the stripe Fermi surface has been 
observed both in systems which show a striped ground 
state and in systems where long-range stripe order has 
not been detectedi^ As in the heavy Fermion case for 
relatively small z we expect that rik is dominated by the 
incoherent spectral weight. According to our scenario rik 
should resemble the Fermi surface of stripes in mean-field 
regardless of whether stripe order is static or fluctuat- 
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ing at low frequencies. Indeed LDA computations in the 
presence of stripes with magnetic and charge long-range 
order reproduce this Fermi surface.— 

The experimental stripe-like Fermi surface is not flat 
as could be expected for a perfect one dimensional band 
structure but shows some wavy features, which depend 
on the details of residual hopping processes perpendicular 
to the stripes. Remarkably the wavy features are well 
reproduced by the LDA computation both with regard to 
amplitude and periodicity, giving credibility to the idea 
that the static LDA computation provides a snapshot of 
the fluctuating order in the disordered phase (cf. also 
Ref . [13 for a more phenomenological approach) . 

The main point is that the formation of QPs due to 
some coherence effect is a small perturbation for the over- 
all distribution of spectral weight. Thus a careful study 
of 7ik can give precious information on the proximity to 
some ordered phase. By the same token the presence 
of QPs and a Fermi-liquid-likc Fermi surface are not in- 
compatible with two-particle responses (neutrons, opti- 
cal) which show strong features of ordering (like stripes) 
at frequencies above loq- This explains why computations 
of fluctuations on top of stripe phases with long range 
orde r 11 ' 12 explain well optical conductivity and neutron 
scattering data of systems without long-range order. It 
also explains how nodal quasiparticles can coexist with 
fluctuating stripe order. 

In the following section we present a toy model of the 
Kampf-Schricffer-typsi 3 - where the dichotomy between a 
Fermi-liquid-like Fermi surface and a momentum depen- 
dent incoherent part, reflecting the fluctuating order, can 
be illustrated. Although we study self-energy corrections 
self-consistently, the lack of vertex corrections makes our 
computations reliable only in weak coupling where QP 
weights are close to one. Within this limitation one can 
show that the described scenario holds. 

In Sec. |TT] we present the model, while in Sec. IHII we 
present some numerical results. A discussion of the re- 
sults and our conclusions are reported in Sec. HVl 



II. THE MODEL 

In order to substantiate the above ideas we consider a 
system of electrons coupled to a dynamical order parame- 
ter which can describe charge ordering (CO) fluctuations 
or spin ordering (SO) fluctuations or both. To fix ideas 
we consider CO fluctuations described by an effective ac- 
tion 

S = -g 2 }2 dTl dT 2Xq(n ~ T2)p q (ri)p_q(r 2 ). 
q Jo Jo 

(3) 

In order to simplify the calculations wc consider a 
Kampf-Schricffer-type model susceptibility 13 which is 
factorized into an u>- and q-dependent part, i.e. 

Xq M = W(ioj)J(q) (4) 



The (real-frequency)-dependent part is W(cu) — 

), with F UOt r being a normal- 
ized lorentzian distribution function centered around ujq 
with half-width T, i^ ,rM ~ T/[(u - lo ) 2 + T 2 ]. The 
momentum-dependent part in D dimensions reads 

M is a suitable normalization factor introduced to keep 
the total scattering strength constant while varying 7 oc 
£ _1 (with £ the CO/SO correlation length). To simplify 
the treatment and to make the effect of fluctuations as 
clear as possible, we will mostly consider the case of spa- 
tially coherent fluctuations. Although formally the cor- 
relation length £ diverges, the ordering is not static and 
we will show that this is enough for the spectral function 
to converge to the Fermi-liquid Fermi surface at zero fre- 
quency. 

The infinite correlation length case is described in D 
dimensions by, 

1 ° 

n=i 

In previous works the spatially-smeared version of Jjj (q) , 
Eq. ([5]), was considered to describe the kink in 
the electron dispersions^ and the (still experimentally 
controversia l 15 ' 16 ) isotopic dependence of these disper- 
sions^ The susceptibility x contains the charge-charge 
correlations in the case of charge fluctuations and spin- 
spin correlations in the case of magnetic fluctuations. 

If only the dynamical part W(uj) were present in 
Xq(iw), one would have a bosonic spectrum B(u>) = 
teuah(w/(kT))F u)0t r((d) which is a "smeared" version of 
the Holstein phonon considered in Ref. [H. The crucial 
feature of the susceptibility (j4|) is the substantial momen- 
tum dependence, which describes the (local) order for- 
mation and reflects the proximity to an instability with 
broken translational symmetry. 

The static limit F UOl r{v) = 8{y) together with an infi- 
nite charge-charge correlation length (7 — > 0) as the one 
considered in Eq. ^ reproduces mean-field results for a 
long-range phase. 13 ' 19 The static limit of (d]) has been the 
object of an intense activity based on the idea, pioneered 
in Ref. l20|, that different types of slow order-parameter 
(OP) fluctuations (S(M or COi 9 -) can be treated as 
classical fluctuations with time-independent correlations. 
The equivalence of such static degrees of freedoms with 
quenched impurities has more recently been formalized 
and cast in field-theoretical language^ Although they 
allow for (nearly) exact solutions, the main drawback of 
these static approaches is that they do not allow for the 
aimed separation of energy (i.e. time) scales since to 
justify the static character of the OP fluctuations, they 
assume that their relaxation time top is much longer 
than the inelastic scattering rate of the electrons r e , 
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top 3> r e £L Here we consider instead 
Fv ,r(v) = S(u - uj q ) 



(7) 



representing a dynamical fluctuation oscillating at a fre- 
quency ujq and therefore averaging out on timescales 
larger than tqp ~ 1/wo- 

In the present limit the problem has also a simple 
Hamiltonian formulation which wc introduce for later 
use: 



H 



+9 



& c L c k<r 

kcr 



Q°Q 



b-Q b - 



(8) 



kcr 



[ C k-Qcr c kcr(6q + 6_ Q ) + 4 +Qo .c kff (^ Q + 6 Q ) 



where c ko . creates a free Fermion and fe^g creates a 
bosonic collective mode excitation and £k = £k — M with 
£k the non interacting dispersion relation and the chem- 
ical potential. 

QPs at energies much larger than lu see an essen- 
tially static fluctuation and modify their dispersion as 
in the mean-field calculation mentioned above. For them 
the fluctuations are static and, if the momentum depen- 
dent part J(q) is strongly peaked around the ordering 
wavevector Q, they are scattered like in the presence of 
a long-range symmetry breaking. Their dispersions are 
modified accordingly. On the other hand, for QPs at en- 
ergies in a shell of width u>o around the Fermi level the 
order parameter fluctuations are far from being static 
and have not enough energy to excite them. For these 
QPs the system is still in a uniform Fermi-liquid metal- 
lic state. Therefore, if their dispersion is followed near 
the chemical potential [m, a large FS fully preserving the 
Luttinger volume is found. In the following we explicitly 
calculate the band dispersions and the FS to show this 
physical mechanism at work and the resulting dichotomy 
of the FS. 

Neglecting vertex corrections we iteratively solve at 
zero temperature the following coupled set of equations 
for the self-energy and the Green's function 

E(k,iw) - -|-^Xq(^)G(k-q,^-ip) (9) 



G(k, iu) 



1 



iu — £k — S(k, iu) 



(10) 



To the best of our knowledge our self-consistent treat- 
ment goes beyond previous solutions that remained at 
the perturbative level. The latter can be obtained at 
lowest order by replacing the Green's function in Eq. © 
by the non-interacting G . For our specific model of the 
susceptibility the self-energy becomes: 



E (1) (k,iw) = 



9 



/(&+q) 



i - m 



k+Qj 



k+Q 



w 



£k+Q . 

(11) 

where the superscript (1) denotes the order of the iter- 
ation and /(£k+Q) = Q(— £k+cj) is the occupation num- 
ber. Inserting this self-energy into Eq. (fTU)) leads to 



a Green's function G' 1 ' (k, iu) which displays now two 
poles and which can be used to compute E^ 2 ) (k + Q, iu>) 
(and thus G^ 2 ^(k + Q,icj)) and so on. The detailed pro- 
cedure is described in the Appendix. 

For the sake of simplicity, we first implement our nu- 
merical analysis in one dimension and then we describe 
the case of two dimensions, which is more relevant for 
layered materials. 



III. NUMERICAL RESULTS OF 
FLUCTUATING ORDER 

A. Numerical analysis in one dimension 

Coherent order parameter fluctuations 

We first present the results for a one-dimensional 
model with a commensurate ordering wavevector Q — 7r 
corresponding to CO/SO with a doubling of the unit cell. 
In particular we consider a band of QPs — — cos(k) 
(we choose a unit hopping parameter t — 1 and a uni- 
tary lattice spacing) coupled via a CO/SO mode given 
by Eqs. dU) and ([6j> . Fig. [2] reports the simplest case in 
which only two poles are considered in the Green's func- 
tion given by G^ 1 -* for a generic filling (n — 0.67: here 
and in the following densities are defined as total num- 
ber of particles per site). This is the standard mean-field 




FIG. 2: (Color online) Upper panel: Two-pole band structure 
in one dimension for g = 0.5t, u>o — and particle density 
n = 0.67. The width of the curve is proportional to the 
weight of the state and energies are measured with respect to 
Ef- Lower panel: Momentum distribution curve. 

solution, with a doubling of the unit cell and a folding 
of the bands resulting in a double FS. In the absence 



5 



of interaction the Fermi momentum would be given by 
k F = nn/2 = 1.05. The dispersion relation instead dis- 
plays two Fermi points and a very different Fermi surface 
volume, "violating" the Luttinger theorem. 

The situation is very different for a dynamical fluctu- 
ation. In Fig. [3] we report the spectral function with 
luo = 0.3 and moderately weak coupling g 2 /u)o = 0.5. 
Also in this case shadow bands appear so that the elec- 
tronic structure has bands close to the location of the 
bands in the broken symmetry state. The shadow bands, 
however, do not reach the Fermi level and one has only 
one FS point at k — kp. Therefore the Luttinger theo- 
rem is satisfied since the divergent pole in G(kp,uj = 0) 
from the static solution now turns into a step so that 
G(k,uj = 0) < for k > k F . 

The spectral function can be written in the Lehmann 
representation on the real axis as: 

A(k,o M ) = ^\(^\ct a \4,$)\ 2 5(u;-E^ + E») 

(12) 

A(k, w < M ) = ^|<^- 1 | CkCT |^)| 2 %-< + ^- 1 ) 

V 

A(k, uj) has structure at the energies of the excitation 
of the system with one added particle (w > ^) or one 
removed particle (ui < fi). The result of Fig. [3] can be 
understood from the low energy excitation of the system 
when 3 = 0. These are listed in Table U and the cor- 
responding excitation energies are shown in Fig. [¥] For 
g = all the weight is in the main band labeled £k be- 
cause the matrix elements in Eq. (|12p vanish when one 
boson is present. The effect of a finite g is to give some 
spectral weight to the "shadow" band at £k-Q±^o and to 
introduce some level repulsion when the bands cross. The 
important point is that the shadow band never touches 
the Fermi level but it is separated from it by the energy 
to create the bosonic excitation. We associate the main 
band with the QP band and the shadow band with the 
incoherent spectral weight. Clearly close enough to the 
Fermi level only the QP band exists. 
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Removal uj < 




State 




Momentum 


c_ ko -|0o r ) 




k 


c_ k+QCT fcl Q |^) 


C-k+q — Wo 


k 



TABLE I: Low energy excitations in the g = limit. The 
central column shows the excitation energy. Notice that 

£-k = Ck. 

Notice that the shadow bands are quite similar to the 
case of a symmetry-broken state (but for a shift in energy 



of order ojq = 0.3, which is small on the scale of the 
hopping t = 1). 

Because the QP states at the Fermi level are negligibly 
affected by the scattering with the dynamical ordering 
mode, the FS (here represented by points) is preserved 
and no shadow branches appear at low energy. There- 
fore already this simple weak-coupling case shows that 
an ordered state would be inferred from the presence of 
shadow bands at high energies, while the low-energy QP 
are characteristic of a uniform state. 
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FIG. 3: (Color online) Upper panel: Two-pole band structure 
in one dimension for g 2 /uo ~ 0.5, ojo ~ 0.3 and electron 
density n — 0.57. The width of the curve is proportional 
to the weight of the state and energies are measured with 
respect to Ef. Lower panel: Momentum distribution curve. 
For k > 2 the occupation number has been scaled by a 
factor 5 to enhance the visibility. 
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FIG. 4: (Color online) Schematic electronic structure for g = 
0. 
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In Fig. [3] we also report the momentum distribution 
rik- From this quantity one can see that some weight 
is indeed present at high energy via the shadow band 
appearing at momenta above k* « 2.3. This is a sig- 
nature of some ordering in the system. However, since 
the CO/SO is dynamical the dispersion discontinuously 
jumps from above to below the Fermi level without pro- 
ducing an additional branch to the FS and with no vio- 
lation of the Luttinger theorem. In this last respect we 
explicitly checked that the Green function at zero fre- 
quency is positive for k < hp, it changes sign through a 
pole at kp and stays negative all the way up to n. The 
Fermi momentum kF is the same that one would have 
in a non-interacting tight-binding system with the same 
number of particles kF = mt/2 with n = 0.57 in this 
case. 

The above picture remains valid if additional poles 
are considered in the Green function provided the QP- 
CO/SO-mode coupling is in a weak-to- moderate regime. 
Fig. [5] reports the case of three poles for the same param- 
eters of Fig. [3J The increased allowed number of poles 



incoherent spectral weight retains a strong momentum 
dependence. 

Upon increasing the QP-CO/SO mode coupling g one 
eventually enters a regime where our non-crossing pertur- 
bative scheme neglecting vertex corrections breaks down 
and the Luttinger theorem is violated. This situation 
is reported in Fig. [6l where a second branch of the FS 
appears due to the strong bending down of the band at 
momenta k w k* . In this case the Green function at zero 
frequency G{tu = 0) changes sign twice (and diverges) at 
the two Fermi momenta kpi and kpi while it changes 
sign passing through zero at k = k* , where the jump in 
the dispersion below fi signals a divergence of the self- 
energy (TTj) . One can check explicitly that the volume 
corresponding to a positive G(u> = 0) is larger than the 
one given by occupied states in the non-interacting sys- 
tem thus Luttinger theorem is not satisfied. We believe 
this is an artifact due to the lack of vertex corrections 
which become important as the coupling is increased. 
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FIG. 5: (Color online) Upper panel: Six-pole band structure 
in one dimension for g 2 /ljq = 0.5 and uio = 0.2. The electron 
density is n = 0.55. The width of the curve is proportional 
to the weight of the state and energies are measured with 
respect to Ep. Lower panel: Momentum distribution curve. 
For k > 2 the occupation number has been scaled by a 
factor 5 to enhance the visibility. 

better represents the shift of spectral weight at high en- 
ergies induced by multiple shadow bands. Nevertheless, 
the finite frequency of the mode protects the low-energy 
QPs from being scattered and leaves the low-energy spec- 
trum unaffected: The FS is still formed by just two points 
despite the appearance of several shadow bands at high 
energy typical of a CO/SO state. The effect of the self 
consistency is to "blur" the shadow bands but still the 
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FIG. 6: (Color online) Two-pole band structure in one di- 
mension for g 2 /ljo = 2. and u>o = 0.3. The electron density 
is n = 0.83. The width of the curve is proportional to the 
weight of the state and energies are measured with respect to 
Ef. Lower panel: Momentum distribution curve. 

The apparent success of this simple theory on illus- 
trating the spectral function of a system with fluctuating 
order is encouraging. However, a different problem arises 
if one considers the half-filled case. Without long range 
order we would expect a metallic state if we were able to 
solve the model exactly. On the contrary commensurate 
scattering in the perturbative solution produces the un- 
physical result of a gaped FS in contrast with the fact 
that no true broken symmetry is present in the system. 
On the other hand it is not surprising that our perturba- 
tive approach fails when our singular interaction %(q, cj) 
connects two degenerate states at the FS as it occurs in 
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this case. We believe this failure is due to the lack of 
vertex corrections which we expect to suppress the scat- 
tering at the Fermi level. We have tested this idea phe- 
nomenologically by assuming that the QPs at the Fermi 
level are protected against this singular scattering by a 
momentum dependence of the coupling of the form 



ff(£k,£k4 



<7tanh(— ) tanh( £k+q ' 



(13) 



One technical remark is in order here. Perturbatively one 
could introduce on the r.h.s. of Eq. (flU)) the bare QP dis- 
persion £k- However, at moderate-strong couplings the 
QP dispersions are substantially modified by the cou- 
pling with the modes to dressed QP dispersions e~k. To 
suppress the scattering of the dressed QPs near the true 
FS, one must insert in (|13[) the renormalized dispersions 
£k- This need of a self-consistency scheme considerably 
complicates the calculations and led us to consider only 
simple symmetry breakings (cell doublings). 




FIG. 7: Two-pole band structure in one dimension for 
g 2 /ujq = 0.5 and wo = 0.3 at half-filling n — 1. In (a) the 
bare vertex g is used, while in (b) the phenomenological ver- 
tex g of Eq. (|13|l is used. 



mostly effective when the exchanged momenta upq are 
larger than the typical exchanged energy ujq. This is the 
case here, where the exchanged momenta are peaked at 
sizable Q's. 

Figs. [7(a) and (b) illustrate the effect of this vertex 
correction in the half- filled case n — 1, where it plays 
a crucial role to restore the metallicity of the system. 
For the previously considered cases of generic filling we 
find that the vertex corrections play a minor role and the 
results obtained with g differ little at moderate coupling 
from those reported above. However, for strong coupling 
the phenomenological form of the vertex Eq. (fT3"| again 
prevents the system to violate the Luttinger theorem as 
e.g. in case of the result shown in Fig. [51 



Finite range fluctuations 

Since in realistic systems the dynamical fluctuations 
have a finite coherence length and a finite lifetime, we 
also investigated the case of a one-dimensional QP band 
exposed to fluctuations having finite T and 7 in i^ 0) r 
and in Eq. Obviously the finite extension in space 

and time of the fluctuations produces broadening in the 
spectra. As it is natural, one finds that the broadening 
in the momentum direction is ruled by 7, while T rules 
the broadening along the energy axis. As it clearly ap- 
pears in Fig. [51 the broadening of the spectra does not 
spoil the essential feature of the coupling with a dynam- 
ical mode (cf. Fig. [3] for 7 = T — but same parameters 
otherwise) : The shadow bands persist as high-energy sig- 
natures of a (local) order, while the FS stays unchanged 
and Luttinger theorem is obeyed provided the mode is 
sufficiently narrow in energy (7 < luq) so that one can 
neglect the "leakage" of weight down to the Fermi level. 




Although the form (|13[) is adopted on a purely phe- 
nomenological basis, we like to remind that several mi- 
croscopic calculation o 22 ' 23 show that in strongly corre- 
lated systems the coupling between QPs and phonons is 
severely suppressed at low energies. This suppression is 



FIG. 8: Two-pole band structure in one dimension for 
g 2 /wo — 0.5 and wo = 0.3 at generic filling n = 0.57. Here a 
finite inverse coherence length 7 = 0.02 and time F = 0.1 are 
used. 
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Pairing effects 

We further explore our one-dimensional model to in- 
vestigate the effects of a particle-particle pairing. To this 
purpose we introduce in our bare QP band structure a fi- 
nite gap A, which, as it is usual for this type of pairing, is 
tight to the FS. Accordingly the Bogoliubov particle-hole 
mixing occurs near the FS and the QP band opens a gap 
centered at the Fermi level. The branch below the Fermi 
level bends down giving rise to a maximum occurring at 
a momentum coinciding with the Fermi momentum of 
the unpaired QPs. Fig. [5] displays this effect for our typ- 
ical parameter set at a generic filling. In particular one 
can see that the particle-particle pairing only modifies 
the low-energy states over scales of order A, while leav- 
ing unchanged the high-energy states, which still display 
the clear effect of the dynamical ordering. This find- 
ing shares close resemblance with recent scanning tun- 
nelling experiments in cuprates^, where an energy scale 
Ao (analogous to our loq) separates the low-energy Bo- 
goliubov coherent QPs from the high-energy excitations 
carrying information of a charge ordering tendency. The 
one-dimensional character of our toy model helps in sepa- 
rating these different excitations in momentum space. A 
more detailed analysis of this issue would require mod- 
eling specific features of the cuprates and is beyond the 
scope of the present paper, but we find this preliminary 
one dimensional finding rather suggestive. 




.1 i 1 1 1 1 — 1 1 — i Q 

0.5 I 1.5 2 2.5 3 

k 



FIG. 9: Two-pole band structure in one dimension for 
g 2 /ioo = 0.5 and ujq = 0.3 at generic filling n = 0.4. Here 
a finite inverse coherence length 7 = 0.02 and time F = 0.1 
are used as well as a finite pairing gap A = 0.1. 



B. Numerical analysis in two dimensions 

Once the main effects of the dynamical mode exchange 
have been presented in one dimension, we illustrate the 
two-dimensional case. For the bare electron dispersion 



we use: 

£k = — 2t(cos(k x ) + cos(ky)) — At 1 cos(k x ) cos(k y ) 

For concreteness the parameters have been chosen to re- 
produce the FS of LSCO, taking t = 342 meV, t'/t = 
—0.2, and /1 = — O.lt. For the mode frequency and cou- 
pling, we choose u>q = 50 meV and g 2 = 1.5o>o, while we 
consider only coherent fluctuations with vanishing broad- 
enings (r — 7 = 0). To avoid too many band foldings, 
which could complicate the analysis, we choose a order 
given by a simple cell doubling in the x direction. This 
gives rise to Q = (±7r, 0). Of course this does not corre- 
spond to the charge/spin modulation observed in LSCO 
cuprates, but it better fits our simplicity and illustrative 
purposes (on the other hand it corresponds to the spin 
ordering in FeAs stoichiometric compounds^-^). 

Figs. [10] (a) and (b) report the two-dimensional FS 
obtained by integrating the spectral density in Eq. ([1]) 
over a small (W; = 10 meV) and a large (Wh = 100 
meV) energy range respectively. Only two poles have 
been considered in the recursive Green's function and 
the QP-mode coupling is dressed via the 'self-consistent' 
vertex corrections [i.e. corrected by the interaction it- 
self, see the remark after Eq. (fl~3|) ]. These corrections 
are required here because at the typical fillings we con- 
sider, the two-dimensional FS has branches connected by 
the critical wavevectors (the so-called "hot spots" in the 
framework of superconducting cuprates). The states con- 
nected by Q would display a gap opening similar to the 
(spurious) one obtained in one dimension at n = 1. For 
this reason we decided to phenomenologically suppress 
this exceedingly strong effect which we attribute to the 
simple perturbative treatment of the singular interactions 
in the model. 

In the case of integration over the low-energy states 
only, one obtains a FS closely resembling the one of un- 
perturbed free QPs'. On the contrary, upon integrating 
over a broad energy range, the "FS" appears folded and 
closely tracks the one expected for a system with a static 
long-ranged broken symmetry. Therefore, also in two di- 
mensions we see that coupling the QPs with dynamical 
modes preserves the QPs from a strong rearrangement of 
all the states (including those near the FS). Notice that 
the vertex corrections are relevant in this regard only 
around the hot spots, but do not prevent the generic for- 
mation of shadow bands or related features appearing in 
n k . 

IV. DISCUSSION AND CONCLUSIONS 

The results of the previous sections clearly indicate 
that dynamical CO/SO fluctuations, at least in a not 
too strong coupling regime, do not produce substantial 
effects on QPs around the FS. Therefore the absence of 
low-energy signatures of CO/SO is compatible with a dy- 
namical type of order. Formally we have used a Kampf- 
Schrieffer-type 13 approach extended by a self-consistent 
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0.5 I 1.5 2 2.5 3 



FIG. 10: Two dimensional FS obtained from the momentum 
distribution function rik in Eq. Q). The spectral function has 
been integrated over: (a) a low-energy range Wi = 10 meV; 
(b) a high-energy range Wi = 100 meV. 



evaluation of self-energy corrections. However, due to 
the lack of vertex corrections our computations are re- 
stricted to weak coupling therefore QP weights are close 
to one and the incoherent weights are small. To have 
the reverse ratios of weight we need to go to strong cou- 
pling, which is unfortunately not feasible. Therefore we 
can only speculate on how the spectral function will look 
like. 

We expect that as the coupling is increased some fea- 
tures of our computations will persist. For example the 
fact that the momentum dependent incoherent part of 
the spectral function is separated by a bosonic excitation 
from the Fermi level is expected to be quite robust in 
strong coupling. Indeed one can construct approximated 



excitations as in Table U but with heavily dressed QPs 
rather than with particles which suggest that a similar 
physics will be at play. More complicated excitations 
may give spectral weight closer to the Fermi level but 
for those the momentum dependence will be completely 
washed out providing a structureless background. 

The fact that the incoherent part should resemble the 
electronic structure of the ordered state follows from con- 
tinuity arguments. For example, as one crosses a transi- 
tion where the spin gets disordered at low energies one 
does not expect dramatic changes in the overall distribu- 
tions of weight in the spectral functions. Those will be 
determined to a large extent by the short range correla- 
tions which may change very little across the order dis- 
order transition. In this sense a faint QP peak emerging 
at low energies, which is certainly a dramatic perturba- 
tion from the point of view of the Fermi liquid properties 
becomes a small perturbation from the point of view of 
the spectral weight distribution. Thus we expect this 
resemblance to persist in the strong coupling limit. 

The scenario we propose is of obvious pertinence in 
cuprates, where standard ARPES experiments usually 
report a Fermi-liquid Fermi surface or at least well de- 
fined nodal particles which at first sight are incompatible 
with vertical stripes. On the other hand other experi- 
ments display the signatures of dynamic order like the fa- 
mous hourglass dispersion relation observed in cuprates 26 
even when CO/SO is not detected) 27 ' 28 This hourglass 
dispersion has been explained computing the fluctuations 
on top of an ordered ground state which will not display 
nodal quasiparticlesiii Our analysis reconcile, at least 
qualitatively, this apparent contradiction: fluctuating or- 
der can coexist with a Fermi-liquid like Fermi surface (in- 
cluding nodal states) and still behave at energies larger 
than loq as a stripe ordered state. 

It is also interesting that recent scanning tunneling 
microscopy data— reveal the simultaneous existence of 
low energy (Bogoliubov) QPs and high energy excitations 
which break translational invariance which thus closely 
agrees with the prediction of our analysis. In particular, 
Fig. [5] clearly shows the simultaneous presence of a su- 
perconducting gap with Bogoliubov QPs at low-energy 
and of high-energy spectral features, which clearly carry 
information on a spatial order. For these high-energy 
states the construction of a local density of states natu- 
rally displays the structure of the staggered order related 
to the scattering wavevector Q. 

As we mention in the introduction Zhov* 7 ^ and collab- 
orators find different Fermi surfaces when analyzing the 
spectral function integrated in a narrow window around 
the Fermi level or in a large energy window. Fig. [TU] il- 
lustrates the weak coupling version of these effect with 
the small window Fermi surface being Fermi liquid like 
and the large window pseudo-Fermi surface showing the 
effects of fluctuating order. 

The dichotomy between low and high energy spectral 
features as discussed in the present paper is closely re- 
lated to another dichotomy observed in ARPES experi- 



10 



merits, namely the one between nodal and antinodal QPs 
in underdoped cuprates* 3 ^ For stripe-like charge order pa- 
rameter fluctuations the characteristic scattering vectors 
are of a comparable magnitude than the separation of 
the Fermi surface segments at the antinodal regions. As 
a consequence (and quantitatively analyzed in Ref. 
within a similar formalism) the scattering is much more 
effective for antinodal states than for nodal ones since 
nodal quasiparticles get protected by momentum conser- 
vation in addition to the dynamical protection discussed 
here. Antinodal states will be affected by the vertex cor- 
rections for which we have only a weak coupling guess 
and therefore require further study but it is conceivable 
that spectral weight at the Fermi level will be suppressed 
in these regions due to transfer of spectral weight to the 
shadows. 

The ascription of the dichotomy between nodal and 
antinodal QPs to low energy char ge o rder parameter fluc- 
tuations is also discussed in Ref. [l9| within a frustrated 
phase separation scenario. For charge scattering along 
the Cu-0 bond direction one also obtains a FS where the 
spectral weight is confined to the nodal regions and thus 
is compatible with a momentum dependent pseudogap in 
the sense of an anisotropy in along the FS. 

One difficulty in interpreting data in the cuprates is the 
presence of disorder. It is possible that the system has 
an ordered ground state but, due to quenched disorder or 
a complex energy landscape — , it never becomes ordered. 
As for a structural glass, the system has long range in 
time positional order but a structural factor characteris- 
tic of a disordered state. In this case scattering probes 
will fail to detect order even though the charge is not 
fluctuating and is inhomogeneous. Local probes, how- 
ever, should be able to detect the ordering but those are 
more scarce and more difficult to interpret. Both static 
disorder and our picture will predict similar results at 
high energies but will strongly differ at low energies were 
only within our picture one recovers a uniform Fermi liq- 
uid. It is not clear at the moment which effect prevails 
in different systems. 

There are cases in which CO has been detecte d 32 ' 33 i 34 
and still the stripe Fermi surface computed in LDA^ and 
measured from the does not show up close to the 
Fermi level. For example no visible shadow FS's appears 
in Ref. H, where the FS of a CO LBCO sample only 
displayed Fermi arcs due to a particle-particle pseudo- 
gap. It is possible that in this case the system charge 
orders because the commensuration with the underly- 
ing lattice provides a pinning potential, which helps to 
stabilize the charge but the spin is still quantum disor- 
dered. Indeed spin order breaks a continuous symme- 
try whereas commensurate CO, as usually observed in 
cuprates, only breaks a discrete symmetry. Therefore 
the former is much more fragile than the latter. Thus we 
propose a state in which the charge is ordered but for the 
spin our dynamical picture applies. One can still wonder 
why one does not see shadows due to the CO. One should 
keep in mind that the CO is a minor perturbation to the 



electrons compared to SO. In a SO state the effective po- 
tential seen by the electrons oscillates on a scale of U. 
In order to check the effect of CO alone without SO we 
have computed the FS in the presence of a charge mod- 
ulation similar to the one observed. The result is a large 
Fermi liquid like FS with weak distortions 3 ^. These may 
be related to features that only very recently have been 
resolved for stripe-ordered LSCO codoped with EuM- It 
is only at moderate energies, above loq for the magnetic 
excitations, that the mixed CO-SO fluctuations produce 
shadow features in the spectra and give rise to "crossed 
FS's" like those observed in LSCO in Refs.0S 

A momentum dependence of the so-called "hump" fea- 
tures at relatively high energies (of order 0.1 eV) has also 
been detected in P^S^CaC^Os+a 38 . One can speculate 
that scattering with a magnetic mode related to SO as 
described here is responsible for this effect. (The authors 
offer a related explanation involving the (tt, tt) spin res- 
onance). An analysis of nk similar to the one carried 
out for LSCO in Refs. 0H could also help in discerning 
whether this is just a remnant effect of the proximity 
to the insulating antiferromagnet or due to fluctuating 
stripes. 

In conclusion, in the light of the analysis carried out 
in this paper, we propose that cuprates are affected by 
finite-energy spin/charge fluctuations related to proxim- 
ity to ordered phases. Whether CO is actually realized in 
a static way is rather immaterial from the point of view 
of low-energy ARPES spectra due to the weakness of the 
charge modulations. Only when the energy is above the 
spin and charge fluctuating scale the electronic spectrum 
is sizably affected. At these energies, however, the detec- 
tion of CO-SO-related pseudogaps (not tied to the FS) is 
quite hard due to the largely incoherent character of the 
spectral lines. Only passing through this dynamical 
tendency to order becomes visible. 

We acknowledge interesting discussions with C. Castel- 
lani, C. Di Castro. M. G. acknowledges financial support 
from MIUR Conn 2005 n. 2005022492 and M.G and G. 
S. from the Vigoni foundation. 



APPENDIX A: ITERATIVE SOLUTION OF THE 
MODEL 

In this Appendix we give a detailed description of the 
iterative solution of the two coupled equations (|9ll0p . In- 
serting the lowest-order self-energy EW(fc, ioS) (Eq. (HU) 
into Eq. (jTUJ) leads to a G^'(k, iu) which displays two 
poles and which can be used to compute YS 2 >{k + Q, iu>) 
(and thus G^'(k + Q,iu>)) and so on. This procedure 
therefore generates the series 

G (1) (fc, iuj) -> G (2) (k + Q, ioj) -» G (3) (k, ioS) 

G (4) (k + Q, ioj) -> G (5) (k, iuj) ► G ( "- 1} (jfe + Q, iw) 
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where G^'^ik + Q,uj) has a n-pole structure that can 
be represented as 



r 

G^(k + Q, l u;) = J2 — 



(») 



(k)} 



(Al) 



Consequently the n-th order for the self-energy is given 
by 



„(«) 



iu ± luq — (k) 



and the sign of loq in the denominator depends on 
whether f(Ei n) {k)) = 0,1 (cf. Eq. (fTT|)). 

Examining the equation for the n-order Green's func- 
tion (iu - e k - s(»+i))G(«+i) = l it turns out that the 
solution can be conveniently obtained by solving the ma- 
trix equation: 



[iul - M] £ = I 



(A3) 



(A2) with 



/ s k ga i ( l \k) 
ga[ n) (k) E[ n \k) T^o 
ga^\k) 



M 








\ga ( n\k) 



(A4) 



E^\k) T uoJ 



and we can identify the (n)-order Green's function as the 
(ll)-element of the matrix K.. Denoting by T the trans- 
formation which diagonalizes M, thereby yielding (n+1) 

eigenvalues Ei n+1 \k), the (n)-order Green's function is 
thus obtained as 



n+l 



[Tis 



i no - e ( ; i+1] \k) 



which also yields the new weights a^ 1 (k) = T\ 



(A5) 



This scheme allows for a systematic evaluation of the 



Green's function up to some given order (n). In the first 
step, the matrix M_ is of order 2x2 with (k) = 1 (cf. 

Eq. ©) and E { n\k) = £fe+Q. Diagonalization yields the 
weights and energies of G^\k,iuj) which can be used to 
construct the matrix for (k + Q, iuj) (where the (11)- 
element of M in Eq. (|A4|) is Eu+q) and so on. This 
procedure creates spectral functions with more and more 
poles, however, in case of a not too strong coupling, it 
will converge in the sense that the weight of newly created 
poles becomes smaller and smaller so that the series can 
be cut at the desired accuracy. 
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